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changes drastically when we study "vortex removed" configurations. 
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1. Introduction 

The material of this talk has a substantial overlap with that presented in the talks of N. Cundy, 
C. Gattringer, T. De Grand, J. Greensite, J. Hetrick, I. Horvath, Y. Koma, S. Solbrig, B. Svetitsky, 
S. Syritsyn. 

We study SU(2) lattice gluodynamics and in our calculations we use the massless overlap 
Dirac operator ||l|] : 

D = ^fl + -£=y A = Dw--, (1.1) 

where A is the Wilson Dirac operator with negative mass term. Anti-periodic (periodic) boundary 
conditions in time (space) directions were employed. It turns out that for SU(2) gluodynamics 
the optimal value of p parameter is 1.4. Furthermore, we have used the minmax polynomial ap- 
proximation to compute the sign function sign(A) =A/VaA^ = 75sign(//), where H = 75A 
is hermitian Wilson Dirac operator. In order to improve the accuracy and performance about one 
hundred lowest eigenmodes of H were projected out. Note that the eigenvalues of ( |1.1D lie on the 
circle of radius p centered at (p , 0) in the complex plane. In order to relate them with continuous 
eigenvalues of the Dirac operator the circle was stereographically projected onto the imaginary 
axis [^. The information about statistics and number of gauge field configurations used can be 
found in 

Below we discuss localization properties of eigenmodes of the overlap Dirac operator. A 
natural measure of the localization is provided by the inverse participation ratio (IPR) which is 
defined as follows (for review see for example Ref. [||]). Let 

Pa W = vlixWxix) , £paW = 1 , 

X 

where Yx {^) is an eigenmode of the overlap Dirac operator in a given gauge field background with 
virtuality X, Dyx = ^¥x- Then for any finite volume V the IPR is defined by 

h=VY,plix), (1.2) 

X 

and characterizes the inverse fraction of sites contributing to the support of p^{x). Note that for 
delocalized modes p^{x) = l/V and hence = I, while an extremely localized mode, p^{x) = 
5v,xo> is characterized by = V. Moreover, for eigenmodes localized on a fraction / of sites (so 
that sup p;t =Vf = f V) we have Ix=V /Vf = 1 //. 

2. Dependence of the IPR on the eigenvalue 

In Fig. [1] we show IPR for various lattice spacings at fixed physical volume. Two different 
effects are clearly observed. First, the IPR grows as function of A„ with A„ —>■ 0, for fixed total 
volume and fixed lattice spacing. Second, for sufficiently small, fixed values of A„ the IPR grows 
with the lattice spacing tending to zero. An alternative way of representing the growth of IPR with 
a ^ is provided by the data on Fig. ^ Here we plot the ratio, R of the densities of the eigenmodes 
with IPR> 5 and with IPR< 5 for two values of the lattice spacing a. The data demonstrate that 
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Figure 1: IPRs for low lying eigenmodes at various lattice spacings and fixed physical volume. The "mo- 
bility edge" Xcr ~ 150 - 200 MeV is clearly seen. 
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Figure 2: The ratio of the densities of the eigenmodes with IPR> 5 to that with IPR< 5. 



the observed growth of the IPR with a — > is due to 'typical' eigenmodes, not due to a few modes 
with a huge value of the IPR. 

Another qualitative feature of the data is existence of a kind of a 'merging point', Xmerge ~ 
150 — 200 MeV. Namely, for A„ > Xmerge the values of IPR practically do not depend on the eigen- 
value Xn and the lattice spacing a. Note that existence of a point where all the curves giving the 
dependence of the IPR on A„ is also seen in the data on localization properties of eigenfunctions of 
a test color scalar particle [^]. Our data indicate scaling in physical units, at least an approximate 
one, of the position of the merging point. 

Thus, our data on the localization properties of the low-lying eigenfunctions of the Dirac 
operator indicate existence of a variety of phenomena. If we look for analogies into solid-state 
physics (reviewed in, e.g., [^]) then, probably, one can find an analog to the growth of the IPR 
as function of A„ with diminishing A„. Indeed, it is known that if one approaches the mobility 
edge, Xcr ' from above, then the IPR blows up as (A„ — Acr)^" where a > 0. Note that this 

' Let us remind the reader that the mobility edge is defined in terms of dependence of the IPR on the total volume of 
the system, Vtot- Namely, (IPR) ~ V,o, for V,o, <x below the mobility edge, A„ < X^y and {IPR) ~ const for Viot — > °° 
above the mobility edge, Xn < Xcr- Our data on the dependence of the IPR on the total volume can be found in the 
original paper \m. 
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Figure 3: The histogram of the distribution of the topological charge for full and vortex removed (VR) 
configurations. 

simple dependence holds for a fixed, large total volume Vtot and the singularity is smoothened due 
to finiteness of the total volume. Now, for chiral fermions one expects that the mobility edge is 
Xcr = 0. Then, the growth of the IPR with A„ ^ indicated by our data (see Fig. [I|) is reminiscent 
of existence of the critical exponent a in the non-relativistic case. 

Since the data exhibit strong dependencies of the IPR on various variables, detailed analysis of 
these dependencies goes beyond the scope of the present study (some details, though, can be found 
in the original paper [^). Here, we confine ourselves mostly to discussion of qualitative effects. 
From this point of view dependence on the ultraviolet cut off, that is the lattice spacing a, seems to 
be most exciting. We will discuss this point further in Sect. 4. 

3. Removing of central vortices 

The vortex-removing procedure |0] drastically changes the properties of the vacuum: con- 
finement and chiral condensate disappear. Below we describe how the vortex removing acts on 
topological susceptibility of the vacuum, distribution of the Dirac operator eigenvalues and local- 
ization properties of the eigenmodes. On Fig]3| we show the distribution of the topological charge 
for full and vortex-removed configurations. It occurs that topological susceptibility, < /V >, 
vanishes after removing of the center vortices. The topological charge of the configuration was 
defined by the index of overlap Dirac operator Q = — n- and n± is the number of exact zero 
modes with positive (negative) chirality. 

On Fig.^(left) we show the quantity < Tip (A) >5 which is related to chiral condensate through 
the Banks-Casher relation [||], {^r^f) = —np{Xn 0), where p(A„ 0) is the density of the (de- 
localized) zero modes in the limit of infinite volume V. For finite lattice volumes these modes are 
near-zero and this is the reason of the limit A„ —>■ 0. It is seen that for the vortex removed config- 
urations there exist the gap in the distribution of the near-zero eigenvalues, thus chiral condensate 
vanishes for vortex removed gauge fields. 

On Fig.^ (right) IPR for full and vortex removed configurations is shown. After vortex remov- 
ing all modes (except of a few surviving zero modes) become delocalized, (IPR) ~ 1. 
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Figure 4: < ;rp (A ) > 5 (left) and IPR (right) for full and vortex removed (VR) configurations. The average 
is taken over 200 gauge field configurations at j3 = 2.5 on lattice 16^ • 18. 



Thus, we can conclude that the non-trivial localization properties of the eigenfunctions of the 
Dirac operator are directly related to the central vortices. Note that similar observations were made 
by using other types of lattice fermions, see in particular and references therein. 

4. Discussions 

The most striking feature exhibited by the data is a strong dependence of the IPR on the lattice 
spacing a. The data were fitted by 

h =bo + h-a^-^ , (4.1) 

where bo i,d are constants to be determined from the fit. The IPR data for the low lying modes 
strongly suggest that 

d ^ . (4.2) 

In other words, the volume occupied by the eigenfunctions tends to zero in physical units as a'^. 

Note that the a dependence of the IPR was suggested first by the data of Ref. [|^]. The lattice- 
spacing dependence exhibited by our data is considerably stronger than that reported in Ref. [p^]. 
In particular, in Ref. [^] the values of IPR were changing, say, between {IPR) ~ 2 and (IPR) ~ 3 
as a result of varying the lattice spacing. In our case, the highest values of IPR are approaching the 
value of 20. At this moment, it is difficult to clarify the reason for this difference. First, we use 
different types of fermions. Also, the color group is different, as well as the lattice action used. 

The strong a dependence of the IPR is in sharp contrast with scaling, in physical units, of the 
density of states and of the chiral condensate (defined through the Casher-Banks relation). Thus, 
the emerging picture is that the size of the eigenfunctions shrinks strongly with a — > while the 
eigenvalues remain stable. In other words, the data suggest that the phenomenon of fine tuning is 
quite a common feature of the Yang -Mills dynamics ^. 

^The standard definition of the fine tuning, which goes back to the Higgs physics, is as follows. If a relativistic 
system has a typical size ro then the energy levels of order A ~ 1/ro are 'natural'. If, on the other hand, A <C l/^o the 
level is 'fine tuned'. 
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Note that the phenomenon of the fine tuning was observed first by studying the properties of the 
lattice monopoles and vortices, for areview see, e.g., [11]. However the definition of the monopoles 
and vortices involves projected fields which makes theoretical interpretation more difficult. The fine 
tuning seen in the data on the fermion localization is defined in explicitly gauge invariant terms. 
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